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ABSTRACT 

Partition functions of eigenvalue matrix models possess a number of very different descriptions: as matrix integrals, as solutions 
to linear and non-linear equations, as r-functions of integrable hierarchies and as special-geometry prepotentials, as result of the 
action of W-operators and of various recursions on elementary input data, as gluing of certain elementary building blocks. All 
this explains the central role of such matrix models in modern mathematical physics: they provide the basic "special functions" to 
express the answers and relations between them, and they serve as a dream model of what one should try to achieve in any other 
field. 
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Matrix model theory [1] studies the integral 
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over N x N Hermitian matrices M as a toy-example of quantum field and even string theory. It is spectacular, 
how much one can learn from this seemingly obvious problem. 

What does it mean to study an integral? 



First, we can simply take it. In this particular case the answer is simple: 



^> ■ Z N (g) ~ (2^ 2 /2 (2) 



and does not look very interesting. However it only seems so. As usual, of interest is not the answer itself, but 
its decomposition, implied, by internal structure of our "theory" . And the more we know about these structures, 
the more interesting decompositions we can obtain. In this particular case we could notice that M = UDU\ 
where D = diagja^} matrix, made from eigenvalues of M , and U is a unitary matrix. Then the same integral 
^vq . is decomposed into two - over unitary matrix U and over N eigenvalues Factoring away the volume Vn 

of the unitary group, we obtain: 

1 r N M T/JV r 
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This is already a somewhat non-trivial decomposition, because 

( 27T )N(N+l)/2 

v "-^r (4) 

what is a considerably more complicated expression than the original (2). 

Second, to study an integral in QFT sense means to treat it as measure, and consider all possible correlators. 
This means that of interest is not the (1) itself, but the averages 

/ . ■ \ ftrM* 1 ...trM ik e'^ tTM2 dM 
d. i, = ( tr M 1 . . . tr M tk ) = i , — 5 

or even their connected counterparts, like 

Cf° nn = Ca - QCj (6) 
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This is already a far-less-trivial problem, and looking at the very first examples one immediately observes an 
emergency of new structure: 

C (AO = N 
C 2 (N) = gN 2 

C 4 (N) = g 2 (2N 3 + N)~ 2(gN) 3 + g 2 (gN) 
C 6 (N) = g 3 (5N 4 + 2N 2 ) ~ h{gNf + 2g 2 (gN) 2 

(7) 

The fact that each correlator is a polynomial (not a monomial) in TV is encoded in the idea of loop expansion. 
The fact that all coefficients are integers signals about connection to combinatorics and is encoded in the idea 
of topological theories. 

Third, if we move in the direction of string theory, we need not just correlators: we need generating 
functions. For the set of Ci lt ... t i k there are two obvious options: 
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where V(z) = YlT=o I^r^- One can a ^ so introduce the connected resolvent 

m 

p^ n {z} = 1 [ly{z i )\ a gZ Il (12) 

i 

The fact that correlators Cj where polynomials in TV is now expressed in the genus expansion of the free energy 
and connected resolvents: 



F{t\g,N}^Y,9 2p ^F P {t\gN} (13) 

p=0 



and similarly 
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Already at this stage something highly non-trivial shows up. This becomes clear from a look on the first few 
resolvents: 
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They all are meromorphic (poly) differentials on a Riemann surface 

S : y 2 = z 2 - A{gN) (16) 

which is called the spectral curve. 

According to the string-theory approach, from this point we should move far enough in a number of different 
directions. 

Other phases. As soon as we introduced the generating function Z{t}, we can start treating it non- 
perturbatively. This means that tk are considered not just as infinitesimal expansion parameters, defining a 
germe, but as the coupling constants, and study what happens when they take finite (or even infinite) values. 
Then Z{t} defines a partition function of a family of theories, called non-perturbative partition function. This 
partition function can be re-expanded not only around the Gaussian point, but around any background potential 
V(M) — ^2 k TkM k . Partition function (particular branch of it) then becomes also a function of parameters 
Xfc, which parameterize the moduli of the spectral curve. Phase transitions take place when the genus of the 
curve changes - it is controlled by the number of extrema of the background potential. The study of these 
dependencies is the subject of Seiberg-Witten theory [2], in matrix-model context the corresponding field is 
sometime called the theory multi-cut solutions or of the Dijkgraaf- Vafa phases [3] . The particular branch of 
partition function is also known as CIV prepotential [4] . The most interesting feature of this prepotential are 
Seiberg-Witten special-geometry equations, describing dependence on the moduli by introducing very special 
"flat" coordinates a k instead of Tk- 

a k = § Ak n 

(17) 

aa k JB k 

and the role of the Seiberg-Witten differential on the spectral curve is presumably played by the 1-point resolvent 
Q(z) = p^\z) [5]. The system of interrelated multidensities p^'"^ can m f ac t be built in a universal way for 
arbitrary Seiberg-Witten family of spectral curves £ - this procedure is now known as AMM/EO topological 
recursion [6] and has surprisingly many applications. Whenever partition function can be reconstructed in this 
way, this signals about the matrix-model hidden behind the scene - and there are already numerous examples, 
when recursion works, but the matrix model is not yet found. 

Various limits. Non-perturbative partition function has a huge variety of different limits and critical 
behaviors in the vicinities of all its numerous singularities. The standard large- N, genus-zero and multiscaling 
limits are just the examples. Related problem is the study of convergency properties of various perturbative 
series. All this is very important in applications and constitutes, perhaps, the biggest parts of traditional 
matrix-model theory. 

Other observables. In string-theory paradigm there is no special preference for any obvious choice of 
observables. Instead of the correlators Cj of the monomials trMj one could study those, say, of the "Wilson 
loops" tr(e sM ), and form many other generating functions, different from (8) and (9), like the celebrated 
Harer-Zagier exact 1-point function [7] 

.0 S, C*" 1 ) 11 (1-A)((1-A)-(1 + A)^) 



and Brczin-Hikami integrals 



I„u iSi ,. f, , S A TT (m» - Uj)(ui - uj + Si - Sj) 
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The number of integrals here is k, not N, as in (3). In fact, these two subjects are unexpectedly closely related 
[9]. Harer-Zagier functions capture contributions from all genera - they differ from (8) by a kind of Padc 
transform and allow to put under the control the divergence of perturbative genus expansion. Instead they 
hide all the information related to spectral curves and Seiberg-Witten equations - but are capable to provide a 
closed expression for the Seiberg-Witten differential fl(z) = p^\z). Unfortunately, they are much more difficult 
to study than the resolvents. 
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Alternative formulations. For non-perturbative partition functions integrals (be they matrix or func- 
tional) provide only a description of particular phases: or, in worst case just the perturbative germes at particular 
points. More adequate are formulations in terms of D-modules or r-functions, characterizing partition functions 
as solutions to linear or quadratic equations respectively. It is still unclear, how general is the existence of 
quadratic (integrability-theory) structures and if higher non-linearities can also be relevant. At this moment, 
the "matrix-model T-functions" - usually, KP/Toda- functions, satisfying also a linear string equation, and, as a 
corollary, a whole infinite set of linear "Virasoro constraints" [10] are the most profound special functions, en- 
countered in modern mathematical physics. They are natural for presentation of quantitative results in various 
fields of string theory, and their investigation is one of the primary purposes of modern science. 



Integrability and TU-representation. Emergency of non- linear (integrable) relations, like [11] 

d 2 log Z N _ Z N+1 Z N -x 



(20) 



for (8), is so non-trivial and so universal in string theory, that it can be considered as one of the main features 
of non-perturbative physics - still very mysterious. One should look for adequate ways to characterize these 
structures. Non-trivial T-functions can be made from the "trivial" ones by integrability-preserving transforms, 
described in terms of the VF-operators, which move the points in the Universal Grassmannian, parameterizing the 
space of the KP/Toda (free-fermion) r-functions. In other words, a matrix- model r- function can be considered 
as a result of the "evolution", driven by cut-and-join (W) operators from some simple "initial conditions" [12]: 

Z{t} = e™T {t} (21) 
For (8) this VT-representation looks as follows: 

Z N {t} = e^\^(at a bt b g^ + (a + b + 2)t a+b+2 ^)\e^ (22) 



a, b 



Generalizations. According to string-theory paradigm, one should not just embed original model in a 
set the similar ones by exponentiating all naive observables, one should also deform everything else, including 
the discrete parameters. In application to matrix models this means that starting from (1) one should not 
just switch from quadratic to arbitrary potential, not just treat N as one of parameters, but also substitute 
Hermitian matrices by others: unitary, orthogonal, symplectic, belonging to exceptional and other Lie algebras, 
to generic tensorial categories etc etc. Of all this the most far-going so far are extensions to unitary matrix 
models [13] and to /3-ensembles [14]. In all cases one expects to find all the relevant representations: not only 
through traditional integral formulas, but also as D-modules, as r-functions, through lU-operators, through 
topological recursion which start from peculiar spectral curves, through Harer-Zagier-type recursions. Some 
results in these directions exist, but they are far from being exhaustive. 

External fields and dualities. Another generalization is inclusion of external fields. The simplest possi- 
bility is to switch from (1) to 

Z(y\A) = e-i trA2 f e- ^ m2 +trMA (dot MY dM (23) 

Determinant is introduced here to make the dependence on A non-trivial, and we also changed the notation for 
the size of the integration matrix. This is done on purpose, because if this function is considered as a function 
of the variable p k = tr A~ k , it is actually independent of n. In this way one defines Kontsevich matrix models, 
eq.(23) is the Gaussian one, for properties of generic Kontsevich models see [15]. Really remarkable is the 
duality between (23) and (8): 

Z(N\A) ~ Z N {t k } (24) 

provided t k = \pu = jftr A~ k . In fact, this duality [16] can be used in the derivation of Brezin-Hikami formulas 
(19), which, in turn have non-trivial generalization [17, 9] to at least the cubic Kontsevich model. 

Unification. Duality between Gaussian Hermitian and Kontsevich models is just an example of interrelation 
between two a priori different matrix models. The goal of string theory is to unify in a similar way all quantum 
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field theories, and in particular, this applies to unification of all matrix models. Unification does not mean 
solving - that problem belongs to the field of non- linear algebra [18], which studies formulas like 



// 



dxdy e *x 2 +bxy+dy 2 „ = D ~V^ 

^Aad - 6 2 2,2 



dxdy e ax +ox y+cxy +ay ~ D 2|3 

D 2 | 3 = 27a 2 d 2 - b 2 c 2 - 18abcd + 4ac 3 + Ab 3 d (25) 

(in general ordinary discriminants D N \ r control singularities of integral discriminants). Unification means that 
all seemingly different non-perturbative partition functions cither arc interrelated (by dualities), or are all 
reductions of some larger partition function (arise at particular loci in the extended space of time- variables) , or 
are all composed from some elementary building blocks. It turns out that the last, most promising, possibility 
can be true, at least in the world of the eigenvalue matrix models. Namely, at least all the Dijkgraaf-Vafa 
partition functions can be obtained by a universal gluing procedure from a few basic elements [19]: 

k 

Z{t} = e e ]JZ^{t^} (26) 

i=l 

where U is bilinear in derivatives over tW-variables. This gluing procedure is closely related to AMM/EO 
recursion [6] and can be considered as one of its most profound implications. The role of the elementary building 
blocks ZW play several important matrix models which posses a sphere with punctures as their spectral curves: 
the Gaussian Hermitian model, the cubic Kontsevich model and the Brezin-Gross-Witten model [20] . 



Applications. Matrix model theory has infinitely many applications in all branches of science, far beyond 
pure mathematics, string theory and even physics. Still, it deserves mentioning a few relatively new examples, 
concerning the abstract fields of research, in order to illustrate once again the influence of matrix model intu- 
ition on our understanding of basic problems. These recent applications also emphasize the role of the character 
calculus — one of the most important matrix-model-theory technical methods. Moreover, matrix models them- 
selves are not present very explicitly, what are discovered are the typical structures and relations, pertinent for 
matrix-model partition functions. 

The first subject is Hurwitz theory [21]. Today it is clear that this is basically the story about the algebra 
of cut-and-join operators, which are well known in matrix model theory 

i v ' 

They are labeled by Young diagrams R = {n > r 2 >} and have Schur functions (the GL(oo) characters) xq[M] 
as common eigenfunctions: 

W rX q[M] = <Pq(R)xq[M] (28) 

while eigenvalues ipQ (R) depend on a pair of Young diagrams and are essentially the characters of symmetric 
group ^(oo). The Hurwitz partition functions describe the sums like 

Q Q \ R J 

and possess many properties, typical for matrix-model r-functions, including the deeply hidden Virasoro- 
constraints, as well as numerous non-trivial generalizations, involving non-commutative "open-string" algebra, 
extending the commutative "closed-string" one formed by the Wr. See [22] for details and references. 

The second subject is the AGT conjecture [23]: the celebrated identity between 2d conformal blocks and 
Nekrasov expansions [24] of the LMNS functions [25], describing instanton expansions of 4/5/6d SYM theories. 
This subject brings together conformal field theory, Seiberg-Witten theory, classical and quantum integrable 
systems [26]. At the core of the story is the special "conformal" matrix model [27], which realizes Dotsenko- 
Fateev representation of conformal blocks, and LMNS functions appear from Selberg integrals, arising in the 
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character expansion of the model. In this language the AGT relation reduces to the Hubbard-Stratanovich 
duality [28], this works perfectly at j3 = 1, but generalization to /3-ensembles remains subtle [5]. What is 
extremely important in this story is that the averages of characters are again characters - and matrix models 
with this special property seem to become more and more distinguished in modern applications. 

The third example is the modern theory of knots [29, 30], which studies extended [31] HOMFLY H®{pk \ q} 
[32] and superpolynomials P®{pk \ q,t} [33]. A very interesting matrix model realization here is long known for 
the underlying Chern-Simons theory, but its generalization in the presence of non-trivial knots is so far available 
only for torus knots [34] and for t = q. In general one expects that the model exists, the measure depends on the 
braid realization B, and the HOMFLY polynomial in representation R is an average of the SL(N) character: 

H% = (xr[U]) b (30) 

Like in the case of AGT relation, one expects that with this measure the averages of characters will be again 
simply re-expanded in characters, and such model will be a useful tool to study the character expansions of 
HOMFLY and superpolynomials, which are responsible for the fast progress in the field in recent months. This 
is indeed the case for the torus knot [m, n]: the measure is given by [34] 

N . ? N 

(■- >,„, - n/ *. n ** ^ *» ^ ( - ) <») 

' ' i=\ J i<j 

and < xr[U] >[ m .n] ~ XR j^fc = j , moreover, like with all Selberg-type integrals, this property persists 
for bilinear combinations of characters. 

The forth example, which deserves mentioning is a very similar Chern-Simons type matrix-model represen- 
tation in the very important ABJM theory [35], describing N copies of M2 branes. The only additional 
complication is that cosh factors are also present in denominators of the Vandcrmondc determinants. Despite 
this complication the model was completely solved in [36] at vanishing times, and the required non-trivial 
behavior <~ A 3 / 2 (instead of the usual <~ N 2 ) of the free energy was reproduced in the large-A limit. 

Note that adequate introduction of time variables, suitable for revealing the linear and non-linear relations - 
in the form of Virasoro constraints and KP/Toda integrability respectively - remains a largely unsolved problem 
in all these examples, despite there are already many signals, that these or very similar structures should exist. 
It is one of the primary tasks of matrix-model theory to study and resolve these mysteries. 
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